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Two different approaches, the direct simulation Monte Carlo (l)SMC) method based on molecular gasdy- 
namics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum 
gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. 
The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. 
Results from the two numerical methods are compared with Rothe\s experimental data, in which density and 
rotational temperature variations along the centerline and at various locations inside a low-density nozzle were 
measured hy the electron-beam fluorescence technique. The continuum approach showed good agreement with 
the experimental data as far as density is concerned. The results from the IJSMC method showed good agreement 
with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation 
parameters, such as the gas/surface interaction model, the energy exchange model between rotational and 
translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC 
method. 


Introduction 

T HE mission performance of satellites and spacecraft such 
as on-orbit lifetimes, payloads, and trip times are sig- 
nificantly impacted by low-thrust rocket engines that arc used 
for the control of altitude and trajectory of the vehicles. An- 
other important factor affecting the mission performance is 
the contamination of sensitive instruments and system com- 
ponents of the vehicles in the plume and backflow region of 
the thrusters. Hence, the understanding of the detailed flow 
structure inside low-thrust rocket nozzles is very important 
for the accurate prediction of the thrust and mass flow levels, 
and also for the precise analysis of the plume and backflow. 

For this type of rocket engine, due to the small thrust level, 
nozzle scales are quite small and reservoir pressures are very 
low. Reynolds numbers of the flow in the nozzle are very low 
and rarefaction effects can significantly alter the internal flow 
structure in the vacuum of the space environment. Under 
these conditions, the flow exhibits strong nonequilibrium ef- 
fects, such as slip at the wall, due to rapid expansion into the 
low-density environment. The fluid experiences continuum, 
transition, and free-molecular flow regimes. Consequently, 
conventional continuum gasdynamics that are based on the 
concept of a local equilibrium may not be adequate, and an 
approach based on molecular gasdynamics is required for the 
analysis of the flow. 

Even though there exist many methods in molecular gas- 
dynamics to analyze low-density gas flows, the majority of 
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the works investigating this type of expanding flow employ 
the direct simulation Monte Carlo (DSMC) method of Bird. 1 
Several DSMC investigations that have dealt with this aspect 
include those by Bird." ' Heuser et al./ Campbell/*’ Nelson 
and Doo, 7 Penko et al. / Boyd et al. / and Zelesnik et al/° 
The DSMC' method is a computer simulation technique to 
solve Boltzmann’s equation by modeling a real gas flow using 
a representative set of molecules. However, due to intensive 
computational requirements, current applications of the DSMC 
method are generally limited to near continuum and rarefied 
flows. Continuum methods are usually much more efficient 
than the DSMC method for high and normal density flows. 
Thus, in the analysis of flows that involve both continuum 
and rarefied flow' regimes, it would be reasonable to combine 
both methods. The simplest utilization of both methods is to 
solve rarefied flow' regimes using the DSMC method by ob- 
taining boundary conditions from the solution of continuum 
methods. This, indeed, is the most widely used method in 
analyzing the flows that involve both continuum and rarefied 
flow' regimes. 

In the present study, two different approaches, the DSMC? 
method based on molecular gasdynamics. and a finite volume 
approximation of the Navier-Stokes equations, w hich are based 
on continuum gasdynamics, are employed in the analysis of 
the internal flow structure in a small converging-diverging 
nozzle. The full Navier-Stokes equations are solved by a low'er- 
upper symmetric successive over-relaxation (LU-SSOR) scheme 
developed by Yoon and Shuen M in the RPLUS code. ,: In the 
DSMC? simulation, the boundary conditions for the solution 
domain are obtained from the results of the Navier-Stokes 
solution. Detailed numerical results from the two approaches, 
such as density and temperature inside the nozzle, are com- 
pared with Rothe’s" experimental data. 

Few experimental data are available for this type of low- 
thrust nozzle, and most data deal with gross characteristics of 
nozzle performance such as thrust levels and discharge coef- 
ficients. This type of data does not provide detailed infor- 
mation regarding the internal flow structure. Rothe’s work is 
probably the only one in w hieh detailed low-density flow prop- 
erties have been measured inside a nozzle using the electron- 
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beam fluorescence technique. In addition to Rothe’s work, 
low-density pitot pressure measurements together with nu- 
merical simulations were performed by Penko et al. N and Boyd 
et al /' They compared continuum and DSMC results with pitot 
pressure data at the nozzle exit plane and at various locations 
in the plume region. Comparison between the DSMC' and 
continuum results was also made by Penko and Boyd for the 
flow inside the nozzle. Campbell, s - h Nelson and Doo, 7 and 
Zelesnik et al. ,n have also considered expanding low-density 
flows using the DSMC method, and compared their results 
with experimental data, but in these investigations no com- 
parison was made with experimental data inside the nozzles. 

The present study reports, for the first time, a detailed 
assessment of the DSMC method for low-density internal flows 
by comparison with experimental data. The validity of the 
numerical methods, the DSMC’ and continuum methods, has 
been assessed by comparing simulation results with detailed 
measurements such as density and rotational temperature at 
various locations inside the low-thrust nozzle. In addition, 
special attention is paid to the DSMC’ simulation parameters, 
such as the gas/surface interaction model and the energy ex- 
change model between rotational and translational modes. 

Problem Statement 

For the analysis of low-density nozzle flow , Rothe's 13 ex- 
periment is chosen as a reference problem due to the avail- 
ability of detailed measurements inside the nozzle. Figure l 
illustrates the geometry of the nozzle used in Rothe's exper- 
iment and in the present numerical analysis. The nozzle is 
made of graphite to reduce optical reflections and to minimize 
backscattering and secondary emission of electrons. The sub- 
sonic and supersonic portions of the nozzle are cones having 
half-angles of 30 and 20 deg, respectively, with longitudinal 
radii of curvature at the throat equal to one-half of the throat 
radius. The maximum area ratio at the exit based on the throat 
area is 66. The computational domain for the continuum method 
consists of the nozzle interior (OAFI), and the region extends 
upstream an axial distance of 10 mm from the nozzle throat 
(OK). The shaded region in the lower portion of Fig. 1 in- 
dicates the DSMC simulation domain. The length of the curved 
nozzle contour (IH) is about 0.5 mm. The simulation domain 
extends to an axial distance of 14 mm from the nozzle exit 
plane (AB), to a radial distance of 14 mm from the nozzle 
lip (FD), and to an axial distance of 14 mm from the nozzle 
exit plane into the backflow region (DE). The inflow bound- 
ary is located at the nozzle throat (OI). Supersonic outflow 
boundary conditions are assumed along the outflow' boundary 
(BC, CE, and EG). The test gas is nitrogen with a stagnation 
temperature of T u = 300 K. The flow conditions are listed in 
Table 1. In the table, the throat Reynolds number, Re, = 
2mtTTfi u R n is based on the viscosity at the stagnation chamber 
condition p. tl . Here, the quantity m is the mass flow rate, and 
R, is the throat radius. The Knudsen number is based on the 
throat diameter and the stagnation chamber condition. 



Fig. 1 Geometry of low-thrust nozzle. 


Continuum Method 

In the continuum method, the Navier-Stokes equations are 
solved by the LU-SSOR scheme 11 in the RPLUS code. 1 - The 
code employs an implicit finite volume LU-SSOR scheme to 
solve the full Navier-Stokes equations and the species equa- 
tions in a fully coupled manner. The LU-SSOR scheme em- 
ploys an implicit Newton iteration technique to solve the finite 
volume approximation of the steady-state version of the gov- 
erning equations. Even though the system of equations is 
formulated in a fully implicit and a fully coupled manner, the 
LU-SSOR scheme requires only scalar diagonal inversion for 
the flow equations. This results in a fast convergence rate, 
and the convergence of the Newton iteration method is as- 
sured by the diagonal dominance of the coefficient matrices 
of the LU-SSOR scheme. The code has been applied to var- 
ious rocket nozzles including H,/() : thrusters, ‘ 4 - 15 nuclear 
thermal rockets, lh and low-thrust nozzles. 17 

In the present study, the inflow total enthalpy and pressure 
are assumed to be constant. At the inlet of the nozzle (KJ in 
Fig. 1), the radial velocity is assumed to be zero and the axial 
velocity is obtained by extrapolation from the interior. The 
temperature is obtained from the inflow total enthalpy and 
the velocity. The pressure and density are obtained from the 
isentropic relation and the equation of state, respectively. At 
the exit of the nozzle where the flow is mostly supersonic, all 
dependent variables are extrapolated from the interior. The 
wall is assumed to be adiabatic and the no-slip condition is 
used at the wall. At the axis of symmetry, the radial velocity 
and the radial derivatives of the other dependent variables 
are set to zero. The calculations are made with a 240 x 60 
grid. A very small uniform axial grid is used from the inlet 
to some distance downstream of the throat, after which the 
grid spacing increases. The radial grid size becomes finer near 
the wall. A detailed description of the continuum method 
used in the present study can be found in Ref. 17. 

Figure 2 shows profiles of flow' variables at the nozzle throat 
obtained from the continuum solution that are used for the 
input boundary condition in the DSMC simulations. Here the 
quantities p, U r , U x , 7 , and .S’ are density, radial velocity, axial 
velocit y, tem perature, and the most probable thermal speed, 
S = \ 2R /’, respectively. The subscript, 0, denotes the stag- 
nation chamber condition. 


Table 1 Flow conditions 


Test gas 

N : 

Stagnation temperature, T u 

300 K 

Stagnation pressure, P„ 

474 Pa 

Wall temperature, T w 

3(H) K 

Reynolds number. Re, 

270 

Knudsen number, Kn 

2.3 x 10 ' 



Fig. 2 Profiles of flow variables at the nozzle throat used for input 
to the DSMC simulations. 
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DSMC Method 

The DSMC method is a popular simulation technique for 
low-density flows, and the DSMC code used in the present 
study is based on the same principles as described in Bird, 1 
together with the variable hard sphere (VMS) model 1 * as a 
molecular model, and the no time counter (NTC) method 14 
as a collision-sampling technique. The code has been devel- 
oped at the NASA Lewis Research Center to investigate var- 
ious low-density flows of gas mixtures in arbitrary shaped flow 
domains, with or without chemical reactions.-" 22 Details of 
the code and computation may be found in Ref. 20. 

The VMS exponent to of nitrogen is chosen to be 0.24, with 
the reference molecular diameter of 4.07 x 10 w ni at the 
reference temperature 273 K. 14 Chemical reactions and the 
vibrational mode are assumed to be frozen. For the calculation 
of rotational energy exchange between the colliding mole- 
cules, the Borgnakke-Larsen phenomenological model 23 is 
employed together with the temperature-dependent energy 
exchange probability of Boyd. 21 However, the temperature- 
dependent energy exchange probability is modified to be con- 
sistent with the experimental data for the rotational relaxation 
of nitrogen obtained by various methods and compatible with 
the VMS model. The modification employed in the present 
study is described below . 

In the Larsen-Borgnakke phenomenological model, a frac- 
tion tf) of translational collisions are assumed to be inelastic, 
and the rest of the collisions are considered as elastic. Hence, 
the fraction </> can be interpreted as the average probability 
of rotational energy exchange for translational collisions. The 
rotational collision number Z w is usually defined by 

Z t< (1) 

where the quantity r A , is the rotational relaxation time, and 
r, is the mean collision time for translational relaxation. The 
rotational collision number is the average number of molec- 
ular collisions that are required for rotational-translational 
energy exchange. Therefore, the average probability of ro- 
tational energy exchange for translational collisions can be 
given as the inverse of the rotational collision number: 

<!> - ( 1/Z*) (2) 



Fig. 3 Temperature-dependent rotational energy exchange proba- 
bility for the VHS model. 

the gamma function. The value of the correction factor is in 
a range of 0.996-1.571. depending on the molecular model. 

Then, the temperature-dependent rotational energy ex- 
change probability suitable for the VHS model is given by 

1 . ,, 

1 ^«.vns V’ 

v,ls 1 + % (T*/ry- + 1(77-74) + 21 (777’) 

(A) 

where the quantity T* is the characteristic temperature of the 
intermolecular potential, and (Z w ) is the limiting value. A 
small numerical correction to Parker's expression noted by 
Brau and Jonkmarr* is included in Eq. (6). Figure 3 shows 
various experimental data 24 37 in a temperature range from 
30-3000 K for N-. corrected according to Eq. (3). The solid 
line represents Eq. (6) with the correction factor from Eq. 
(5) for T* - 91.5 K, Z* = 23.5, and = 0.24. The tem- 
perature-dependent rotational energy exchange probability 
employed by Boyd 21 is also shown in the figure. 

Employing the method introduced by Boyd, 24 Eq. (6) may 
be converted into the following instantaneous exchange prob- 
ability suitable for the VHS model: 


As is the case in Boyd, 24 the rotational collision number 
may be obtained from Parker's expression 2S by determining 
the constants to match available experimental data and the 
more rigorous analysis of Lordi and Mates. 2 '' It should be 
noted here that special attention should be paid in the inter- 
pretation of experimental data due to the use of different 
mean collision times for translational relaxation in presenting 
the experimental data. Also, these data must be corrected so 
that they should be compatible with the VHS model 

4,v„s - F.Z r (3) 

where the quantity Z H VMS is the rotational collision number 
suitable for the VHS model, and b, is a correction factor that 
is given by 


F, = <tA,.viis> (4) 

where the quantity r, VM s > s Ihe mean collision time for trans- 
lational relaxation for the VHS model, if r, is given by tt/jl/ 
(4 E), which is one of the most widely used mean collision 
times,’ 7 it can be shown that the correction factor is given by 

l5tt7rF(2 - w) 

f ' ■ »l <4 - (5 > 

where the quantity a is the correction factor to the first ap- 
proximation of the coefficient of viscosity, 14 and F denotes 


^VHS (Z H ) b ( 1 + 


I (£ +2 to) / A / * \ 1 it 3| 2 
E(£ + 3/2 - <o) \ E ( / 2 


+ 


\'(( + 2 - to) / / 7T 

\'U + 1 - *>) \ E, / \ 4 " 


(?) 


where the quantity ( is the average internal DOF of the col- 
liding molecules, k is the Boltzmann constant, and E ( is the 
collision energy. 

A diffusely reflecting wall with 10% thermal accommoda- 
tion is assumed for the interaction between the gas molecules 
and the wall. For comparison, several gas/surface interaction 
models are also considered, including diffuse reflection with 
full thermal accommodation, specular reflection, adiabatic 
wall, and a combination of the three models. The energy 
exchange between rotational and translational modes when 
the molecules collide with the surface involves quite complex 
physical processes. Also, these energy exchange mechanisms 
largely depend on the kind of gas molecules and the wall 
materials. One of the interesting examples is the interaction 
between a graphite surface and NO or Br 2 . lK In this case, the 
rotational temperature of the reflecting gases never exceeds 
250 K, independent of the surface temperature above 300 K 
due to the strong rotational cooling effect. Due to the lack 
of experimental data between graphite and nitrogen, the wall 
is assumed to be a monatomic gas with an infinite mass for 
the rotational-translational energy exchange of the reflecting 
molecules, and the probability is assumed to be given by Eq. 
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(7). For comparison, full thermal accommodation of rota- 
tional energy at the wall is also considered. 

Results and Discussion 

To present the general idea regarding the differences be- 
tween the continuum and DSMC methods and the overall 
structure of the flowfield, density and Mach number contours 
obtained by the two methods will be considered first. Figure 
4 shows density contours obtained by the two methods. The 
DSMC solution occupies the upper portion of the figure, and 
the continuum solution is shown in the lower portion. The 
density is normalized by the stagnation chamber density p ( „ 
and plotted in a logarithmic scale with base 10. It can be seen 
that the density profiles near the throat exhibit a density ridge 
where the density is higher than on the axis. The density 
humps indicate a weak compression wave originating near the 
throat where the wall curvature changes abruptly. The DSMC 
method predicts a weaker compression than the continuum 
method. In the entire flowfield, the density obtained by the 
DSMC method is lower than that obtained by the continuum 
method. The fluid experiences about two orders of magnitude 
decrease in density across the nozzle, and about 5 orders of 
magnitude in the backflow region at the plane parallel to the 
exit. In some parts of the backflow region, a severe scattering 
of data in the results of the DSMC method appeared if the 
normalized density is lower than 10 ", which made it difficult 
to interpret the results there. Hence, the density contours are 
plotted up to 10 ( \ Some wiggles in the contours are due to 
the scattering of the data, the interpolating schemes used in 
the plotting software, and the larger cell size in this region. 

Figure 5 shows the Mach contours obtained by the two 
methods. Again, the DSMC solution occupies the upper por- 
tion of the figure, and the continuum solution is shown in the 
lower portion. Along the axis, the DSMC method predicts 
higher Mach numbers than the continuum method. This means 
the DSMC method predicts a slightly faster expansion of the 
flow, i.e., lower density and temperature, than does the con- 
tinuum method. It also can be seen that the DSMC method 



Fig. 4 Comparison of density contours, log, 0 (p/p o ), obtained by the 
continuum (lower! and DSMC (upper) method. 



Fig. 5 Comparison of Mach number contours obtained by the con 
tinuum (lower) and DSMC (upper) method. 


predicts the turning of the flow toward the nozzle lip and the 
termination of the sonic line of the internal boundary layer. 
This is consistent with previous observations — experimen- 
tally, w analytically, 4 " and from DSMC simulation results 2 — 
that if the flow is highly underexpanded the effect of the very 
high-pressure gradient at the nozzle lip is to accelerate the 
subsonic portion of the boundary layer. This results in the 
intersection of the sonic line with the nozzle lip. Another 
difference between the two methods is that the DSMC method 
predicts a much thinner subsonic layer near the nozzle wall 
due to velocity slip. The accurate analysis of the flow structure 
in the boundary layer near the nozzle lip is very important in 
the prediction of plume backflow. This is because the flow in 
the nozzle wall boundary layer near the nozzle lip is the origin 
of the plume backflow, and is, thus, the most important factor 
affecting the flow structure in the plume backflow. 

Consideration is now given to the detailed flow structure 
inside the nozzle. In Fig. 6 density profiles along the centerline 
of the nozzle obtained by the continuum and the DSMC meth- 
ods together with the experimental data of Rothe are shown. 
Densities are normalized by the stagnation value p M , and the 
axial distance, which is measured from the nozzle throat, is 
normalized by the throat radius. The results from both the 
continuum and the DSMC methods show good agreement 
with Rothe s experimental data. The maximum uncertainty 
limit of the experiment for the density data along the cen- 
terline reported by Rothe 41 is ± 10%. Figure 7 shows tem- 
perature profiles along the centerline of the nozzle obtained 
by the continuum and the DSMC methods together with the 
experimental data. Temperatures are normalized by the stag- 
nation value T {] . In the continuum method that is based on 
the concept of local equilibrium, the flow is represented by 
only one equilibrium temperature, while the DSMC method 
gives both translational and rotational temperatures. Due to 
the fast expansion, the rotational temperature is alw ays higher 
than the translational temperature along the axis, and it can 



Fig. 6 Comparison of density variation along the nozzle centerline. 



Fig. 7 Comparison of temperature variation along the nozzle cen- 
terline. 
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Fig. 8 Comparison of temperature variation along the nozzle cen- 
terline for P 0 = 209 Pa. 

be seen that the rotational temperature obtained by the DSMC 
method matches Rothes experimental data. It should be noted 
here that the temperature measured by Rothe is the rotational 
temperature. To compare the effect of the rotational energy 
exchange probability, the rotational temperature obtained by 
using the instantaneous energy exchange probability em- 
ployed by Boyd 9 24 is also shown in the figure. It can be seen 
that the introduction of the correction factor, which is con- 
sistent with experimental data and compatible with the VHS 
model, results in a higher rotational temperature and gives 
better agreement with Rothe's experimental data. The effect 
of the correction factor is more accentuated as the flow be- 
comes more rarefied. Figure 8 shows temperature profiles 
along the centerline of the nozzle obtained by the continuum 
and the DSMC methods together with the experimental data 
for the case of P {) — 209 Pa, which is about 2.3 times lower 
stagnation chamber pressure than that of the reference prob- 
lem. In this lower pressure case, a quite different flow struc- 
ture from the higher pressure case can be observed. In the 
higher pressure case, the axial temperatures decrease mono- 
tonically from the throat to the nozzle exit, whereas in the 
lower pressure case the temperatures pass through a minimum 
and then increase toward the nozzle exit due to rarefaction 
effects. That is, as the flow becomes more rarefied, the effect 
of molecule-surface collisions increases while that of mole- 
cule-molecule collisions decreases. From a continuum point 
of view, this means there is an increase of the thermalization 
of the flow energy by viscous dissipation. The rotational tem- 
perature obtained from the present DSMC calculation shows 
good agreement with the measured rotational temperature 
except for the minimum point. All the numerical methods 
predict the location of the minimum at xiR, s 8, while that 
in the experiment occurs at x/R, = 6. The experimental error 
bars in Figs. 7 and 8 are random errors only, which are based 
on the signal-to-noise ratio in the rotational spectra. 13 Sig- 
nificant differences can be clearly seen from the comparison 
with the measured rotational temperature. Although not shown, 
there is no appreciable difference in the translational tem- 
perature due to the change in the rotational energy exchange 
probability. 

Figure 9, in which contours of the two temperatures ob- 
tained by the DSMC method are shown, demonstrates the 
degree of thermal nonequilibrium between the translational 
and rotational modes. The translational temperature occupies 
the upper portion of the figure, and the rotational temperature 
is shown in the lower portion. Temperatures are normalized 
by the stagnation value T {) . The translational temperature 
drops faster than the rotational temperature along the axis. 
A remarkable difference between the two temperatures oc- 
curs around the nozzle lip where the flow rapidly expands 
around the corner. In this region, the rotational temperature 
becomes almost frozen due to the rapid expansion into vac- 
uum. 



Fig. 9 Comparison of temperature contours TIT 0 obtained by the 
DSMC method. Rotational (lower) and translational (upper). 
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Fig. 10 Comparison of density profiles at the plane near the exit. 



Fig. 1 1 Effect of background pressure on the density profile at the 
plane near the exit. 

Comparison between the continuum and DSMC methods 
is now made in the plane near the nozzle exit. Figure 10 shows 
density profiles at xfR, ~ 18.7, 3 mm inside from the exit, 
obtained by the continuum and the DSMC methods together 
with the experimental data of Rothe . Densities are normalized 
by the value at the axis p, , and the r distance is normalized 
by the nozzle radius at the location R sy . It can be seen that 
as the flow comes close to the nozzle exit, the inviscid core 
almost vanishes and a boundary layer occupies the whole cross 
section of the nozzle. The results from both the continuum 
and the DSMC methods show relatively good agreement with 
Rothe’s experimental data. The slight underprediction of the 
density near the nozzle wall in the DSMC result is due to the 
effect of ambient pressure. The effect of the background pres- 
sure is demonstrated in Fig. 1 1 where the DSMC result with 
a finite background pressure is shown. To investigate the 
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effect of the ambient pressure, the pressure along the outflow 
boundaries (CE and EG in Fig. 1) is assumed to be 1/310 P (l , 
which is the actual ambient pressure in Rothe’s experiment. 
Although the boundaries used in the present calculation may 
not be far enough from the exit and the nozzle lip for an 
accurate result, it is believed that this will give a rough esti- 
mation of the effect of the ambient pressure. It can be seen 
that the effect of applying this ambient pressure increases the 
density near the nozzle lip and makes it coincide with the 
experimental data, while it has a negligible effect on the flow 
around the axis. 

In Fig. 12 density profiles obtained by the two methods at 
various locations inside the nozzle together with the experi- 
mental data are compared. The density is normalized by the 
stagnation value. Rothe’s experimental data is calculated from 
his centerline density and cross-sectional densities normalized 



Fig. 12 Comparison of density profiles at various locations inside the 
noz/le. 



Fig. 13 Comparison of temperature profiles at the plane near the 
exit. 



Fig. 14 Comparison of temperature profiles. 



Fig. 15 Effect of wall interaction models on the rotational temper- 
ature profiles at the plane near the exit. 

by the centerline density. As discussed earlier, the density 
profiles obtained by the two numerical methods show humps 
at xlR, = 3.7, while the experimental data is flat. Except for 
the region near the wall and the large hump in the continuum 
result near the throat, the agreement between the calculations 
and the experimental data can be considered fair. The poor 
performance of the numerical results near the wall is a subject 
of future studies. However, it should be noted here that, in 
order to convert beam intensity to gas density, a knowledge 
of the electron beam current at the point of measurement is 
required. 41 This requirement poses a problem near the wall 
due to backscattering and secondary emission of electrons at 
the nozzle surface, and can lower the accuracy of the result 
near the wall. At higher pressure levels collisional quenching 
effects cannot be neglected, and this also can lower the ac- 
curacy of the experimental data near the throat. 

Figures 13 and 14 show temperature profiles at xlR, - 18.7, 
3 mm inside from the nozzle exit, and at xlR , = 13.7, re- 
spectively, obtained by the continuum and the DSMC meth- 
ods together with the experimental data. Temperatures are 
normalized by the stagnation value T u . The stagnation value 
is almost recovered in the translational temperature obtained 
by the DSMC method near the wall, while in the continuum 
method the recovery of the temperature is about 90%. Once 
again, it can be seen that the rotational temperature obtained 
by the DSMC method shows good agreement with Rothes 
measured rotational temperature. In Fig. 15 is shown the 
effect of wall interaction models on the rotational temperature 
at xiR t = 18.7, 3 mm inside from the nozzle exit, in the DSMC 
method. Four different gas/surface interaction models are 
considered: 1) diffuse reflection with full thermal accommo- 
dation, 2) diffuse reflection with no thermal accommodation 
(adiabatic wall), 3) diffuse reflection with 10% thermal ac- 
commodation, and 4) 70% diffuse — 30% specular reflection. 
The DSMC result with diffuse reflection with 10% thermal 
accommodation at the wall shows the best agreement with 
the measured rotational temperature. Although not shown, 
the gas/surface interaction models considered here have a 
negligible effect on the density profile, except for the specular 
reflection that results in a faster expansion of the flow, i.e., 
higher axial velocity, lower density and temperature inside 
the nozzle. 

Conclusions 

A low-density nozzle flow is analyzed by two different ap- 
proaches, the continuum method based on a finite volume 
approximation of the Navier-Stokes equations, which are based 
on continuum gas dynamics, and the DSMC method based 
on molecular gasdynamics. The results are compared with 
Rothe's experimental density and rotational temperature val- 
ues that were measured at various locations inside the nozzle. 
Comparison of results from the two methods with the exper- 
imental data show that the continuum method can provide 
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relatively good results inside the nozzle as far as density is 
concerned- The results from the DSMC method show good 
agreement with the experimental data, both in the density 
and the rotational temperature. It is also shown that the sim- 
ulation parameters, such as the gas/surface interaction model, 
the energy exchange model between rotational and transla- 
tional modes, and the viscosity-temperature exponent, have 
substantial effects on the results of the DSMC method, and 
a proper choice of these parameters is very important for more 
accurate results. These parameters are largely dependent on 
the kind of propellant and surface materials, and more ac- 
cumulation of experimental data bases are required on this 
subject. 
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